Ab-initio description of hole localization and Zhang-Rice singlets in one-dimensional doped 

cuprates 
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We present the first ab-initio band-theory-based description of spin-compensated polarons (known as Zhang- 
Rice singlets) in a hole-doped cuprate, specifically one-dimensional Ca2+ a; Y2- : ECu50io. Zhang-Rice singlets 
are many-particle configurations relevant to the exotic behavior of hole-doped cuprates, stemming from spon- 
taneous charge localization. They appear in our case-study material above a threshold doping, successively 
turning the insulating undoped antiferromagnet into a gap insulator, a singlet-rich metallic paramagnet, and 
finally, a singlet-saturated insulating diamagnet. 
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An outstanding challenge in condensed matter theory is the 
thorough understanding of the electronic structure of hole- 
doped cuprates. A key ingredient of the physical properties of 
many such materials is thought to be a multiparticle electronic 
state known as spin-compensated polaron or Zhang-Rice sin- 
glet (ZRS in the following), first introduced by Zhang and 
Rice llJ] in the context of underdoped cuprate superconduc- 
tors. A ZRS is in essence a localized, positively-charged spin- 
singlet state appearing when (a fractions of) a doping hole lo- 
calizes on one or more O's first-neighbor of the Cu 2+ site of a 
CuC>2 unit (embedded in e.g. a stripe or layer), coupled anti- 
ferromagnetically to the polarized native Cu hole. 

There is ample evidence of the presence of ZRS in doped 
cuprates. Beside their observation [2] in superconductors 
and their relevant role therein yD, experiments J3, H H 01 
indicate that the prototypical one-dimensional (ID) cuprate 
Ca2+a;Y2_2;Cu50io (henceforth CaYCuO) must be highly 
populated by ZRS in a specific doping range. In this paper 
we study CaYCuO by a fully first-principles band theory (the 
pseudo-self-interaction correction approach [8], already ap- 
plied to other correlated cuprates |9J]) theoretically interpret- 
ing the experimental observation of ZRS at low doping, and 
predicting other ZRS related effects at doping levels not ex- 
perimentally attained so far. 

We point out that this is the first instance of a first-principles 
density-functional-theory-based band approach predicting or 
describing the formation of ZRS, which are usually stud- 
ied within Hubbard-like models [ 1 ] with simplified electronic 
structure. While often (perhaps questionably) dismissed as 
inessential in 2D cuprates, a realistic material-specific de- 
scription is unquestionably necessary in the context of ID 
cuprates. By and large, first-principles calculations have long 
been absent from this arena due to their difficulties in treating 
strong correlated cuprates II 1 Qfl . 

ID cuprates are chain-like aggregations of weakly coupled 
chains of CuC<2 units. Compared to other ID systems (e.g. 
GeCu0 3 or Sr 2 Cu0 3 ) CaYCuO has the advantage 01 of be- 
ing hole-dopable in a wide range. Magnetization and elec- 
trical conductivity measurements J3, H @, 01 in the range 
x/n= [0,0.4] (here x is the hole concentration and n the num- 
ber of CUO2 units per chain, and x/n the hole fraction per 



Cu) show that CaYCuO is an antiferromagnetic (AF) wide- 
gap Mott insulator at x/n=0 and a Curie- Weiss AF insulator 
with one-dimensional hopping-conductivity at x/n<0.3. At 
x/n^0.3 the AF susceptibility spreads out in a broad maxi- 
mum, the signature of a magnetic chain populated by a high 
density of spinless sites (ZRS). The conductivity is thermally 
activated with a gap of 0.08 eV. Thus, perhaps surprisingly, the 
system is effectively insulating through all the experimental 
doping range. No data are reported for concentrations higher 
than x/n=0.4. Here we deliver a microscopic description of 
CaYCuO up to full doping x/n=l, explaining the mechanism 
of the observed phase transitions, and uncovering further, not 
yet experimentally accessible doping effects. 

Structure - CaYCuO (FigQ]) is made of ferromagnetic (FM) 
Cu0 2 chains (Cu-Cu separation 3.4 A, Cu-O-Cu angle ~90°) 
running along y and parallel to the (x, z) plane ||4j, |5L Inter- 
calated Ca/Y chains, parallel to the CUO2 ones, act as electron 
reservoir for the CUO2 units. In the synthesized compound 
Ca2+ x Y2_3;Cu50io the ratio of Ca/Y and Q1O2 chain steps 
is 4/5, with n=5. Contiguous chains are weakly coupled anti- 
ferromagnetically (AF) along z, with Tn^28 K at x=0 (16 K 
at x/n=0.2). Coupling along x is negligible. To gain computa- 
tional flexibility, we treat n as a free parameter (as electronic 
properties depend on x/n but not on n) and perform calculation 
for n= 2,3,4. We use a plane-wave ultrasoft-pseudopotential 
method, with 30 Ryd cut-off and 10 x 10 x lOk-point grids; for 
density of states (DOS) calculations, the linear-tetrahedron in- 
terpolation is used on the same mesh. Full details will appear 
in a forthcoming paper. 

Undoped CaYCuO -The calculated orbital-resolved DOS in 
Figf2]show that undoped (x=0) CaYCuO is a Mott insulator 
with a gap of 3.44 eV and Cu magnetic moment of 0.71 hb 
(near-gap bands are shown and discussed in FigHJbelow). No 
experimental gap value is available, but previous studies for 
low-dimensional cuprates [9] suggest our prediction to be sen- 
sible. We obtain the observed magnetic ordering as energeti- 
cally stable, with FM intra-chain exchange-interaction J y =8. 5 
meV and AF inter-chain J 2 =-0.2 meV (~7 meV and ~1 meV 
experimentally). 

Our calculation attests that CaYCuO is a proper Mott- 
Hubbard insulator, i.e. valence band top (VBT) and conduc- 
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FIG. 1: Structure of CaYCuO. FM y-parallel Cu chains couple an- 

tiferromasmeticallv in the z direction. Dashed lines delimit the or- 
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FIG. 2: (Color online) Orbital-resolved DOS for undoped CaYCuO. 
For Cu the contribution of spin-polarized d yz orbital is represented 
by the (red) filled line, while the (black) solid line sums up the other 
four d orbitals. For O, solid and filled lines represent p y and p z states, 
respectively. No Y or Ca DOS appear in this energy range. 



tion band bottom have the same orbital character. All Cu are 
in the 2+ (d 9 ) state with the highest d yz orbital spin-polarized. 
The flat DOS peak located at VBT and separated from the 
lower-lying valence manifold is due to Cu d yz -0 (p y , p z ) 
hybrid bands originating from intra-chain FM coupling. The 
magnetization is about 70% Cu d yz and 15% from each adja- 
cent oxygen. The gap opens between up- and down-polarized 
channels of these p-d hybridized bands. It is worth noting 
that a standard density-functional calculation predicts a non- 
magnetic metal. 

Doped CaYCuO: Fig J5] summarizes the calculated proper- 
ties vs. doping. Fig|4] reports the corresponding band struc- 
tures. We find three markedly different regions correspond- 
ing to hole concentration ranges 0<x/n<l/4, l/3<x/n<l/2, 
and l/2<x/n<l. The x/n=l case stands out separately, as dis- 
cussed below. 

a) Underdoped region - At low concentration (0<x/n< 1/4) 
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FIG. 3: (Color online) Calculated properties as a function of doping. 
Black circles are relative FM chain magnetization per Cu02 unit, 
green triangles are Cu magnetic moment, red squares are the aver- 
aged number of ZRS per Cu02 unit. Going from x=0 to x=l system 
changes from AF Mott insulator (x=0) to AF hopping insulator (HI), 
ZR gap insulator (GI), ZR metal (M), ZR Mott insulator (at x=l). 



injected holes empty progressively the p-d hybridized VBT 
peak, without major changes in the band shape: the injected 
charge disperses along the chain, and no ZRS forms. The 
chain magnetization decreases almost linearly with x, the AF 
alignment remaining unaffected. As for conduction, Fig0] 
shows that at x/n=l/4 (one hole every four CuC>2 units) only 
the very flat (d yz -p y , p z ) spin-majority bands are cut by Ep, 
and then only in the k^, k z -parallel segments M-R and M'- 
M, orthogonal to the chains. Therefore the system is a low- 
mobility hopping-conductive insulator. This region is char- 
acterized (Figj3) by zero ZRS concentration (red squares), a 
smooth linear decrease of relative chain magnetization (black 
circles), and a weak decrease of Cu magnetic moment (green 
triangles). 

b) ZRS gap insulating region - As doping nears the thresh- 
old x/n=l/3, ZRS populate the system: spin-polarized holes 
localized on the oxygens are coupled with the spin-antialigned 
holes on the adjacent Cu, so the total magnetic moment of the 
corresponding G1O2 unit is S=0. This occurs because hole- 
hole Coulomb repulsion exceeds the Cu-O charge transfer en- 
ergy, favoring hole localization. This doping regime exhibits a 
gapped behavior, with transport properties dominated by op- 
tical absorption conductivity. Again we note that a standard 
local-density calculation would not obtain hole localization, 
hence no ZRS. 

FigO shows DOS and magnetization isosurfaces for 
x/n=l/3. At this concentration the holes are fully localized, 
and form one ZRS each three CUO2 steps. Two band features 
signal the presence of ZRS. The first is the flat states of the 
doping-injected holes, strictly localized on oxygens and with 
marked p 2 character, at ^0.4 eV above the VBT. The second 
signature is the ~ 1 e V downshift (from the undoped value) of 
part of the Cu d yz spin-minority bands. An antialigned pair 
of spin-majority p z and spin-minority d yz adjacent holes cou- 
ple to form one ZRS, whereas the d yz spin-minority Cu states 
not coupled to any localized hole (i.e. not engaged in a ZRS) 
remain in the undoped-like position ^3.5 eV above VBT. 
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FIG. 4: Band energies at various dopings: system is hopping- 
conductive up to x/n=l/4, gap-insulating until x/n=l/2, metallic 
above x/n=l/2, gap-insulating at both end-points x/n=0 and x/n=l. 




FIG. 5: (Color online) Upper panel: orbital resolved DOS of CaY- 
CuO at x/n=l/3 (non-ZRS, left; ZRS, right). Holes localize on O 
p z states ~0.4 eV above VBT, and couple to adjacent Cu d yz states 
(at ~2.5 eV). The p z and d yz orbitals (enclosed by ovals) are anti- 
aligned, i.e. the two-particle system is a singlet (ZRS). Lower panel: 
hole density isosurfaces (h=0.027 bohr~ 3 ) for states within 4 eV 
above VBT. White and gray (yellow) surfaces are up and down po- 
larized holes. In blow-up, the ZRS-forming anti-aligned p z and d yz 
orbitals are highlighted by a dashed circle. 



The ZRS is clearly visible in the spin-polarized hole charge 
isosurfaces of Fig(5] Furthermore, ZRS formation shows up 
in the calculated magnetic moments in Table Q] Non-ZRS Cu 
and O atoms keep their undoped-state magnetic moment. In 
the units hosting the localized holes, Cu and O's compensate 
each other quite exactly, forming a ZRS. 



TABLE I: Orbital contributions to magnetic moments. "Cu, O" refer 
to non-ZRS units, "Cu*, O*" to ZRS units. In the latter, the opposite 
magnetizations on d yz and p z states dominate: the corresponding 
Cu02 unit is demagnetized (S~0). 
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As discussed e.g. in J3, El S 0L the low doping-to-ZRS 
phase transition is associated to a change from AF to some dis- 
ordered paramagnetic (cluster spin-glass or spin-liquid) mix- 
ture of ZRS and small residual FM chain segments. Indeed, 
consider in Figj3]the difference between Q1O2 (black circles) 
and Cu (green triangles) contributions to the average magne- 
tization. The former drops linearly through the whole doping 
range, since holes deplete spin-majority states (either O p or 
Cu d) at any x/n. Cu magnetization, instead, decreases lin- 
early only in the low-doped region, but as soon as ZRS appear, 
it remains roughly constant up to x/n= 1 . This means that as x 
increases, ZRS formation demagnetizes more and more CuC»2 
units, progressively destroying the chain FM ordering, while 
leaving Cu magnetic moments at about their undoped values. 
(SQUID susceptibility data [7] also suggest ordered ZRS pat- 
terns, but calculations for these structures exceed our present 
computing capabilities.) 

c) High-doping ZRS-rich metallic region - As x increases, 
the number of ZRS increases as long as the hole concentration 
is large enough to allow the "condensation" of more holes on 
the oxygens. However, as the average hole-hole distance is 
reduced, the charge associated to each ZRS spreads out: at 
x/n=l/2 (i.e. one ZRS per two CuC>2 unit), injected holes no 
longer localize on a single oxygen, but spill out on adjacent 
CuC<2 units. The corresponding states span a ~l-eV energy 
range, much larger than the optical gap of the ZR-GI regime 
and also larger than the undoped VBT DOS peak. As a conse- 
quence, the high-doped (x/n>l/2) ZRS-rich region is weakly 
metallic, although large resistivity should still be expected 
due to the small band dispersion. Since metallic ZRS cap- 
ture a larger hole charge fraction than insulating ZRS present 
at x/n=l/3, at x/n=l/2 some of the chains remain ZRS-free, 
thus causing the ZRS concentration drop visibly in Figj3]in 
correspondence of the metal-insulating transition. 

According to the hole isosurface plots in Figj6] in this 
regime the systems has no recognizable magnetic order, and 
can be seen as a disordered mixture of impure ZRS (i.e. with 
S not exactly zero: due to metallicity O and Cu polarized 
states do not compensate exactly) and CUO2 units with vary- 
ing magnetization. (These features will be illustrated below in 
an overview of the evolution of magnetic ordering and band 
energies vs x/n.) 
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FIG. 6: Hole density isosurfaces (same as Fig.5) at various dopings: 
up to x/n=l/4 the original AF ordering still rules; at x/n=l/3 ZRS 
form, and one each 3 CuC>2 units form a singlet; intra-chain FM 
ordering is lost, and system is cluster spin glass or paramagnetic. 
Above x/n=l/2 more than 50% of CuC>2 has S=0, i.e. system be- 
comes more and more non-magnetic. 



d) ZRS-saturated diamagnet - The above picture holds for 
increasing doping, until at x/n=l (i.e. one hole for each CuC»2 
unit) our calculations predict an exotic ZRS-saturated insu- 
lating regime, where each Q1O2 hosts one ZRS. The system 
is now a collection of singlet states, hence a non-magnetic (in 
fact diamagnetic) insulator, despite Cu magnetic moments still 
close to their undoped values. It would be exciting to push the 
experimental doping limit up to x/n=l and verify the validity 
of this prediction. 

Overview- In Figj6] several frames of hole density isosur- 
faces are shown as a function of doping. This juxtaposition 
highlights the evolution from low-doping intra-chain FM or- 
dering to ZRS ordering: for x/n<l/3 all O and Cu holes are 
parallel-oriented (thus CuC»2 magnetization is larger than the 
Cu magnetic moment) and hole isosurfaces can hardly be dis- 
tinguished from their undoped counterpart at x/n=0. With 
ZRS formation, holes collapse onto O-centered p 2 -shaped 
charges, antiparallel (different color in the Figure) to the 
dj, z -shaped Cu holes. On the other hand, the small hole 
density residing on non-ZRS oxygens remains spin-parallel 
(equally colored) to Cu holes. The change from ZRS insu- 
lating to metallic regime (x/n^l/2) shows up in the irregu- 
lar hole-charge distributions on the oxygens: now even non- 
ZRS oxygens can give intra-chain AF contributions to the to- 
tal chain magnetization (because of ZRS broadening to the 
nearby CuC>2 units). 

The calculated band structure vs x/n in Fig. [4] con- 



firms again the two regime changes mentioned earlier. The 
transition from low-doping hopping-conductive to ZRS gap- 
insulating at x/n=l/3 is characterized by the flat oxygen bands 
at -0.4 eV above the VBT Thus, in the region l/3<x/n<l/2 
we expect transport properties to be dominated by optical ab- 
sorption through this energy gap. Above x/n=l/2 the optical 
gap closes due to ZRS spreading out over the chains. At x/n=l 
the spectrum is gapped once again: one spin-polarized hole 
per CuC>2 is exactly what is needed to empty out the flat p-d 
bands present at the VBT of the undoped system (see Fig|2]l. 

In conclusion, we described ZRS formation and properties 
in a doped ID cuprate using a first-principles band theory. 
We interpreted in detail the observed simultaneous change 
of magnetic and transport properties of CaYCuO upon dop- 
ing, and extended the analysis to dopings not (yet) reachable 
experimentally. The fact that a first-principles band theory 
proved able to describe ZRS formation in a doped Mott insu- 
lator for the first time is an instrumental step for a realistic de- 
scription of a wide class of low-dimensional doped cuprates. 
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